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We investigate the effect of a unidirectional qnenched random field on the anisotropic quantnm 
spin-1/2 XY model, which magnetizes spontaneously in the absence of the random field. We adopt 
mean-field approach to show that spontaneous magnetization persists even in the presence of this 
random field but the magnitude of magnetization gets suppressed due to disorder, and the system 
magnetizes in the directions parallel and transverse to the random field. Our results are obtained by 
analytical calculations within perturbative framework and by numerical simulations. Interestingly, 
we show that it is possible to enhance a component of the magnetization in presence of the disorder 
field provided we apply an additional constant field in the XY plane. Moreover, we derive generalized 
expressions for the critical temperature and the scalings of the magnetization near the critical point 
for the XY spin system with arbitrary fixed quantum spin angular momentum. 

PACS numbers: 


I. INTRODUCTION 

Disorder is ubiquitous in solid state materials. The 
effectd of disorder are often non-intuitive, and hence 
have generated a lot of attention in condensed-matter 
physicsi^— . Disorder can currently also be engineered 
in a controlled way via ultracold atoms trapped in op¬ 
tical lattices subjected to an additional, e.g. optical 
speckle potential^. Disordered systems are often en¬ 
dowed with nontrivial properties, dramatically different 
from those of their homogeneous counterparts. The novel 
quantum phases^i^ and unique phenomena, such as An¬ 
derson localization^, dynamical many-body localization^ 
and presence/absence of thermalization^, and high 
superconductivity!^ are some of the prominent exam¬ 
ples. Quenched disorder in type II superconductors have 
been investigated by using Ginzburg-Landau theory to 
uncover “glassy” properties in the systemic. In partic¬ 
ular, considerable efforts have been dedicated to under¬ 
standing the effects of disorder in spin models, both clas¬ 
sical and quantumi^^— . 

In classical systems with continuous symmetry, it 
has been shown that a random field with the symme¬ 
try of the system may cause significant changes in its 
propertie a^^d^ . A small random magnetic field of this 
kind can destroy magnetization in a classical spin system 
at any temperature, including zero-temperature. E.g., 
two-dimensional XY and Heisenberg models do not mag¬ 
netize in presence of random fields with SO{2) and SO{3) 
symmetries respectivelyi^di (see alsoi^ for the analogous 
effect in systems with discrete symmetry). The effect 
prevails in higher dimensions as long as the random field 
exhibits the corresponding symmetry!^. However, in ab¬ 
sence of appropriate symmetry, these systems exhibit 
spontaneous magnetizatio n^'^d^ . 


Despite the difficulty in dealing with spontaneous mag¬ 
netization and other system characteristics in quantum 
disordered systems, many important results have been 
obtained. A quantum Hall nematic phase has been pre¬ 
dicted in a zero temperature two-dimensional electron 
system that is unstable to weak disorde r^^d° . Collective 
properties of magnetic impurities on a topological surface 
were studied both theoretically^i and experimentally^. 
Experimental and theoretical investigations in solid state 
systems have revealed that alloy disorder can reduce the 
Curie temperature in the system^^. Arbitrarily weak in¬ 
terparticle interactions were shown to destabilize the sur¬ 
face states of topological superconductors in the presence 
of non-magnetic disorder—. The critical behavior and ef¬ 
fective exponents in ferromagnetic quantum phase tran¬ 
sitions of disordered systems were derived in ReL— . 

Recently, a mean-field classical spin model with SO{n) 
symmetry was considered and its spontaneous magneti¬ 
zation was investigated in presence of unidirectional ran¬ 
dom fields!^. The natural question is how a symmetry 
breaking random field affects these systems in the quan¬ 
tum limit. As already mentioned before, this question 
is particularly relevant due to the current accessibility 
of disordered quantum spin models in experiments^. If 
we restrict ourselves to the XY model with a random 
field, and to one dimension, then such quantum sys¬ 
tem of a moderately large size can be investigated using 
the Jordan-Wigner technique^. However, higher dimen¬ 
sional quantum spin models remain intractable due to 
the lack of analytical and numerical techniques, even for 
ordered systems. We work within the mean-field (ME) 
approximation which is often effective in capturing sys¬ 
tem’s properties qualitatively. Numerical schemes like 
exact diagonalization are usually inefficient for even mod¬ 
erately large systems, due to the exponential growth of 
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the dimension of the system’s Hilbert space. Metasta¬ 
bility effects and slow relaxation rates, usually present 
in disordered systems, make other numerical simula¬ 
tion techniques such as density-matrix renormalization 
group, Monte-Carlo approach difficult to apply, particu¬ 
larly in higher-dimensional lattice systems with higher¬ 
dimensional spins. The mean-field approach essentially 
liberates one from these immediate challenges allowing 
for a detailed analysis, that enables us to answer some 
of the key questions. Of course, the price to pay is that 
the mean field approach is not expected to describe the 
details of the critical behavior precisely, except in high 
dimensions. 

In this work, we consider the quantum spin-1/2 
XY model with anisotropic interaction in presence of 
a unidirectional quenched random field. The purely 
isotropic case, i.e., the quantum XY model with van¬ 
ishing anisotropy parameter and vanishing disorder, ex¬ 
hibits a spontaneous magnetization which has circular 
symmetry. The anisotropy breaks the continuous sym¬ 
metry even for the pure system. The pure spin-1/2 XY 
system magnetizes below a certain critical temperature. 
This system still magnetizes when a random field is in¬ 
troduced at a critical temperature which is higher than 
in the system without disorder. We show by means of 
numerical as well as perturbative analysis that the sys¬ 
tem now magnetizes in specific directions, which is either 
along the parallel or the perpendicular directions to the 
random field. The critical temperature in both cases de¬ 
creases with the increase of the random field strength. 
We find that the critical temperature to magnetize in 
transverse and parallel directions show opposite behavior 
with respect to the anisotropy in the system. Specifically, 
with the increase of anisotropy, the critical temperature 
corresponding to the transverse magnetization increases, 
while the opposite happens for parallel magetization. It 
is important to mention that for vanishing anisotropy 
parameter, the continuous symmetry of the pure XY 
system is broken by the introduction of an arbitrarily 
small random field. We also present general expressions 
of the scalings of critical temperature of magnetization 
for the quantum XY spin systems with arbitrary half¬ 
integer and integer spins. 

In addition, adding a constant magnetic field along 
with the random field, we find that the component of the 
magnetization perpendicular to the random field gets en¬ 
hanced due to the disorder, an effect known as “random- 
field-induced order” or “order from disorder”, which has 
also been reported for several other models^i^l. 

The rest of the paper is arranged as follows. In Sec.llll 
we introduce the spin-1/2 quantum XY model in pres¬ 
ence of the random field, and subsequently derive critical 
scaling of the magnetization via perturbative approach. 
We also discuss numerical results obtained within the MF 
approximation. In Sec. uni we demonstrate, both numer¬ 
ically and analytically, the order from disorder phenom¬ 
ena under the influence of an additional constant field. 
In Sec. El we derive the generalized expressions for the 


case of arbitrary integer and half-integer spins. Finally, 
we conclude in Sec. 


II. QUANTUM SPIN-1/2 XY MODEL IN A 

RANDOM FIELD 

We consider the quantum spin-1/2 XY model in a ran¬ 
dom field. Our aim is to study the effect of the random 
field on the magnetization as a function of temperature, 
and to find the scaling of the magnetization around the 
critical temperature. In the following subsection, we in¬ 
troduce the system and the mean-field approximation. 

A. The system and its mean-fleld treatment 

The Hamiltonian of the ferromagnetic quantum XY 
model is given by 

HxY — ^int T (1) 

where 

N 

Hm = - + Jycrl^)]. (2) 

K-il=i 

The coupling constants, Ja, are assumed to be positive. 
They can be further expressed in terms of an anisotropy 
parameter, 7 , as Jx = T(1 -I- 7 ) and Jy = J(1 — 7 ). The 
indices, i and j denote the sites of an arbitrary d dimen¬ 
sional lattice and af = x^y are the Pauli matrices on 
the site. The part of the Hamiltonian in Eq. o due 
to an inhomogeneous magnetic field, Hext, equals 

Hext — C ^ ^ hi • (Ti^ (3) 

i 

where e (> 0 ), a dimensionless parameter that quantifies 
the strength of the randomness, is typically chosen to be 
small. The unidirectional random field, hi, is chosen to 
be hi = rji ■ ty, where rji are independent and identically 
distributed quenched Gaussian random variables with 
zero mean and unit variance, and Cy is the unit vector in 
the y-direction. Within the mean-field limit, as we show 
below, the pure systems governed by the Hamiltonian 
Hint magnetizes even at low-dimensions. This does not 
contradict the Mermin-Wagner-Hohenberg theorem^^, 
which predicts no spontaneous magnetization in one- 
and two-dimensions since the predictions made by the 
mean-field approximation only becomes accurate in 
higher dimensional^ where Mermin-Wagner-Hohenberg 
theorem is not valid. Interestingly, it has been shown 
that a uniaxial random field may help the system to 
magnetize even at two-dimensio n^^i^° . Note that had 
the random field, hi, been chosen to be invariant under 
rotations, the system would not magnetize at any 
nonzero temperature in any dimension d < 413,17,18 ^ 
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FIG. 1: (Color online) System magnetizes in directions par¬ 
allel and transverse to the disorder field. Zero contour lines 
of the and Fy^{m) in Eqs. (1101) [solid-red] and (lllll 

[dotted-blue] for 7 = 0.3, e/J = 0.1 and J13 = 2, respectively, 
as functions of rrix and niy. All quantities are dimensionless. 


Within the mean-field approximation, each spin is re¬ 
garded to be reacting to an average field due to all 
the other spins in the system. Assuming N to be 
the total number of spins in the system, the effec¬ 
tive interaction, replacing the nearest neighbor interac¬ 
tion in Hint, for large N, equals approximately Hint = 




N 




where TOq, = J' is the coupling constant 

and v is the coordination number which depends on the 
geometry of the lattice. Note that mean-field approxi¬ 
mation provides close to the exact description for large 
dimensional lattice systems^. Within the mean-field ap¬ 
proximation, the Hamiltonian, Hxy, is given by 


H = - J{1 + 'y)mxcrx - J(1 - 'y)myCjy - erf - a, (4) 

where the operators TOq are replaced by their average val¬ 
ues denoted by the same symbol, in the canonical equilib¬ 
rium state at absolute temperature T and J = J'v. Now 
in order to monitor the behavior of the magnetization as 
a function of temperature, one needs to calculate the ex¬ 
pectation value of the spin operators, CTq, a = x,y. So in 
the mean-field regime, the magnetization of the system 
governed by the Hamiltonian, H (see Eq. is given 
by 


rria = Av; 


Tr [tTg exp(-/3i?)] 
Tr [exp(—/?iJ)] 


( 5 ) 


where /3 = l/{kBT), where ks is the Boltzmann con¬ 
stant, and Av^[.] denotes the average over the realiza¬ 
tions of randomness. From Eqs. o and m, we obtain a 
coupled set of the following two equations: 


= mx = Avy 


J(1 


tanh(/3fc'^) 


( 6 ) 


and 


= rriy = Avr, 


J{1- j)my -b £7? 


tanh(/3A;'^ 


( 7 ) 


where k^ = ^/(/(l -I- '))mxY + {ey + J(1 - 'y)my)'^. 

Note that the subscripts _L and || classify two distinct 
cases, as would be apparent later. The superscripts e 
and 2 keep track of the strength of the disorder and of 
the value of the spin (1/2) respectively. 


B. Critical point and scaling of magnetization near 
criticality 

The magnetization, m, can be obtained by finding the 
common zeros of the following two functions obtained 
from Eqs. ([5]) and O) 

= Avy[ '^^^ tanh(/3fc'')] - rrio,, (8) 


and 


tanh(/3fc'')] - rriy, (9) 


where k'^ = 1 /(J(1 -|- j)mx)'^ + {erj + J{1 — 'y)myY. Let 
us set rriy; = mcos^i, rriy = msin/)i, and fh = (rriy;, rriy). 

By performing perturbative analysis, we can study the 
magnetization for small e. Taylor series expansion of Eqs. 
([8|) and in e around e = 0 gives 

F^’^im) = Ci + h,e^ + 0(e'^), ( 10 ) 

and 

Fy'^{m) = C2 + ]^b2e^ + 0(e‘^), (11) 

where 

Cl = rrixi—l + tanh(/3fc)), (12) 

k 

C 2 = my(—l + tanh[/?fc]), (13) 

rC 
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bi = 


-3J^mxmll3{l - 7)^(1 + 7 ) 


+7) 


3 J^mxTOy(l — 7)^(1 + 7 ) tanh[/3fc] 


k'^ 


cosh[/3fc]2 k"^ 

Jmx{^ + 7 ) tanh[/3fc] 
k^ 


k^ 


cosh[/3fc]2 

( 2 J^ma;TO^/3^(l - 7)^(1 + 7 ) tanh[;9fc]) 




cosh[/3fc]^ 


(14) 


-3J^myl3{l - 7 )^ 1 3Jmy/3(l - 7 ) 1 3J^my{l - 7 )^ tanh[/3fc] 

k'^ cosh[/3fc]^ cosh[/3fc]^ fc® 

3 Jmy(l - 7 ) tanh[/3A:] 2 - 7 )^ tanh[/3fc] 1 

/c^ k^ cosh[/3fc]^ ’ 


(15) 


Here k = 7 ) 771 ^;)^ + (J(1 - l)myY. 

A contour analysis at this point becomes helpful to 
characterize the behavior of the system, in particular, 
in finding the directions in which the system magnetizes 
(see Fig. [T]). This amounts to identification of the zero- 
contour lines corresponding to Eqs. (fTOjl and (fTTll . The 
intersection points of the zero-contour lines are possible 
solutions of the magnetization. For any given set of pa¬ 
rameters, one immediately finds that the roots of the 
F'’^(m) and Fy’‘^{m) exist only at (j)i = Q or 7r/2. This 
implies that the magetization is either transverse to the 
random field (case 1) or parallel to the random field (case 
11). Note that for e = 0 and 7 = 0, the zero contour lines 
for both the equations would lie on top of each other due 
to the circular symmetry in the system. However, an ar¬ 
bitrary small random field is enough to break this sym¬ 
metry. It follows from the contour analysis that above a 
certain temperature, the critical temperature, the zero- 
contour lines corresponding to Eqs. m and (HU inter¬ 
sect only if rux = ruy = 0, which is a trivial solution. 

In order to find the critical temperature and the scal¬ 
ing of magnetization near criticality, we perform another 
round of Taylor expansions in Eqs. (nnD and m around 
m = 0 to obtain 

E^’^(m) = --^(3 -I- J/3(l -I- 7)(-3 -|- /3^e^))TOCOS(()i 

-r- J^/3^(l -|- 7)^(5 — 8/1^-I- 4/3^e^ cos 2 (j)i)m^ cos^i 

3! 5 

-I-C>(to^), 

(16) 

and 

Fy’‘^{m) = (—1 -b J/3(l — 7)(1 — e^))m sin (j)i 

— J^/3^(l — 7)^(5 — 16/3^ + 4/3^e^ cos2(/)i)m^ sin(/)i 

3! 5 

-bO(m^). 

(17) 

The contour analysis implies that the allowed values of (pi 
are tt/2 and 0. Eor transverse magnetization, i.e., for the 


case 1 with pi = 0, Fy'‘^{m) vanishes (see Eq. (IT71) 1. The 
nontrivial solutions, which solely appear from Eq. (HU), 
are given by 


rnP^ = ±-\/5< 


' 3 +J^(l + 7)(e"/32-3) 
(—5 -b + 7)3 ’ 


(18) 


The critical point is obtained by setting = 0 in 

Eq. (fT^ . We get 


P 


e ,2 

c,_L 


1 

■/(1 + 7) 


3J3(l-b7)3' 


(19) 


Here P is associated with the subscript T, following the 
similar convention in magnetization. 

The magnetization values corresponding to case 11 are 
obtained by setting = 7 r /2 in Eqs. (fTBl) and (fT71) . In 
this case, the function in Eq. HU vanishes. The nontriv¬ 
ial solutions of Eq. HU are given by 





1 +J/3(l-7)(e2/32-l) 
J3^3(l_^)3(4g2y32_l)' 


( 20 ) 


Subjecting Eq. dini) to the constraint wjj’^ = 0 , we obtain 
the following expression for the critical temperature: 

^‘=’11 ^ J(l-7) ^ J3(1-7)3- 

Eollowing the set of Eqs. (ITOl) - (fTO]) l and (l^ - (l^ l. one 
can immediately infer that the effect of disorder is more 
conspicuous if the system chooses to magnetize along the 
direction parallel to the random field, as compared to 
the other possibiliy with a transverse magnetization. In¬ 
terestingly, the findings are consistent with the pictures 
drawn within a classical limiti^ for the isotropic case i.e., 
7 = 0. However, quantitative comparison shows that the 
analysis with classical spins overestimates the effect of 
disorder on the critical scaling. 

Note that one can immediately deduce the scaling ex¬ 
pressions for the magnetizations and critical temperature 
of the isotropic ordered systems by putting 7 = e = 0 
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FIG. 2: (color online) Numerical and analytical results ex¬ 
hibit persistence of spontaneous magnetization in specific di¬ 
rections even after insertion of disorder. Numerical results 
for the magnetization as a function of J/3, in the directions a) 
transverse and b) parallel to the disordered field. Red circles 
correspond to the roots of Eqs. and © with e/J = 0.1 
and 7 = 0.1. Insets: The blue solid lines correspond to the 
analytic solutions derived for small m given in Eqs. m and 
(uni for the same set of parameters. The red circles are the 
numerical results. We find that the numerical and analyti¬ 
cal results agree in the small m regime. All quantities are 
dimensionless. 

in Eqs. dsni) and (ED) ( or equivalently, in Eqs. m and 
(fT^ ). In this case, the solutions form a circle in the XY 
plane. This can be easily understood by following the 
set of Eqs. © and 0 , which become identical for the 
isotropic ordered systems. The symmetry is broken in 
presence of the random field and then the system prefers 
a specific direction of magnetization. 

Moreover, it is clear from Eqs. (HU and m, that crit¬ 
ical temperature at which system magnetizes increases 
with the coordination number, i.e., with the dimension 
of the system, for any value of 7 . 


C. Numerical results 

In the previous section, we derived the expressions for 
magnetization near the critical point. However, away 
from the critical point, where the perturbative approach 
is no longer valid, one has to rely on a numerical simu¬ 
lation to find the roots of the coupled set of equations, 
given in Eqs. m and ©• We use the classical Monte 
Carlo technique for performing a configurational averag¬ 
ing over 1 ], It requires a few thousands of random realiza¬ 
tions in order to obtain converged values. Our numerical 
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FIG. 3: (color online) Magnetization as a function of J/3 for 
different choices of the anisotropy constant, 7 , in the direc¬ 
tions a) transverse and b) parallel to the random field. Gircles, 
squares, and crosses correspond to the magnetization of the 
system with e/J = 0.1 and 7 = 0.2, 0.4, and 0.6, respectively. 
Insets show the inverse critical temperatures as functions of 
7 for e/J = 0.1. All quantities are dimensionless. 

searches show the presence of two kind of solutions, either 
7 ^ 0, = 0 (i.e., case I) or = 0, 7 ^ 0 

(i.e., case II), which is in accordance with previous discus¬ 
sions in the context of analytical perturbative analysis. 

Fig. UKa) exhibits the results obtained by numerical 
analysis for the transverse magnetization, i.e. case I, 
with vanishing //-component and non-zero a;-component, 
for e/J = 0.1 and 7 = 0.1. When the temperature is 
high enough, the system does not magnetize, similarly to 
the case of an ordered system, i.e. e/J = 0. However, 
below the critical temperature (i.e. if /3 > ,5^’^), the 
system magnetizes in the direction transverse to the ap¬ 
plied random field. We see that the critical temperature 
decreases in the presence of the disorder, i.e. that the 
critical point, shifts towards the right in presence 
of the random field. We also find excellent agreement be¬ 
tween the exact numerical results and the approximate 
analytical expression of the transverse magnetization de¬ 
rived using a perturbative approach given in Eq. (1181) (see 
inset of Fig. [IK a)). 

The numerical results for case H with vanishing y- 
component and non-zero a;-component are shown in 
Fig. |2Kb) for e/J = 0.1 and 7 = 0.1. The features of 
parallel magnetization are qualitatively similar to that of 
transverse magnetization. However, we find that the crit¬ 
ical point, /3^’y, as may be expected by now given the an¬ 
alytical results, shifts towards an even higher value com¬ 
pared to the case of transverse magnetization. A closer 
examination of the Figs. |2Ka)-(b) show that the effect of 
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disorder is more prominent in the parallel magnetization 
than in the tranverse one. This is confirmed by the ex¬ 
pressions derived in the small-m regime (see Eqs. (UHl) & 

m)- 

The behavior of the transverse and the parallel mag¬ 
netizations for a given e and selectively chosen values of 
the anisotropy constant, 7 , is demonstrated in Fig.|31 We 
find that the inverse critical temperature, (3^ j_, decreases 
with increasing 7 for the case when the system magne¬ 
tizes in the direction which is transverse to the applied 
random field. The opposite happens when the system 
magnetizes in the direction which is parallel to the radom 
field. The insets of the Figs.|3Ka-b) show the critical tem¬ 
peratures Pc's as functions of 7 . The trends suggest that 
for highly anisotropic systems, the parallel magnetiza¬ 
tion would occur only at sufficiently low temperatures. 
High anisotropy favors transverse magnetization, i.e. the 
system starts magnetizing in the transverse direction at 
comparetively higher temperatures. 


III. ORDER FROM DISORDER: RANDOM 
FIELD QUANTUM XY MODEL IN THE 

PRESENCE OF AN ADDITIONAL UNIFORM 
FIELD 

Until now, we have seen that the spontaneous magne¬ 
tization in the system persists, albeit only in a restricted 
set of directions, even in the presence of a disordered 
field. Is this still true when there is an additional con¬ 
stant field? In this section, we consider this question and 
show that not only does the spontaneous magnetization 
persist—disorder can now help one of the components 
of the magnetization to achieve an enhanced value com¬ 
pared to the ordered system. 

We first consider the case in which the ordered XY 
model is subject to a constant magnetic field, h. The 
mean-field Hamiltonian, Hh-, governing the system in this 
case is given by 

Hh = - [j ((1 + 'Y)mxcrj: -b (I - 7)771^(7^) + h ■ (22) 

The constant field h lies in the XY plane, 
i.e. h = {hx,hy) = (hcosx, hsince) with magnitude 
h, where 0 < /i < I, and phase x, with —-kI2 < x < -k 12. 
In the presence of the constant field, the mean-field 
equation for the magnetization is obtained replacing 
H by Hh in Eq. (0. The system now has no critical 
temperature, as there is always a unique solution at any 
value of 13. 

Let us now investigate the effect of a random field, erf, 
on the system. The mean-field Hamiltonian Hh.e can be 
written as 

Hh,f. = - {j((1 + l)'mxOx + (1 - l)'m.yay) + h- a + eWy^ 




FIG. 4: (Color online) Order from disorder. Effect of the ran¬ 
dom field in presence of a constant field. The X-component 
of magnetization, rux, as a function of J/3 for (a) 7 = 0.05 
and (b) 7 = 0.1. Red pluses show rrix for the case when the 
XY model is subjected to a constant field h with h/J = 0.3, 
and X = 7 r/ 3 . The blue solid and the green dashed lines rep¬ 
resent nix when the system is treated with random field of 
strength e/J = 0.1 and e/J = 0.2, respectively, along with 
the constant field (the corresponding Hamiltonian is given 
in Eq. m)- The insets show blow-ups of the same for a 
smaller range oi jp. The enhancement of rrix in presence of 
the disorder field uncovers a “random field induced order”. 
Comparing panels (a) and (b), we observe that for fixed h/J 
and X, the order-from-disorder phenomenon appears as the 
anisotropy parameter is cranked down. All other quantities 
are dimensionless. 

where we assume the random field to be directed along 
the y-axis. Replacing again H by Hh.c in Eq. , we 
obtain two coupled equations, which we solve to find fh. 
As may be expected, the solution for the magnetization 
is again unique. 

Our numerical calculations show that the magnetiza¬ 
tion m and the y-component of the magnetization vector 
my are reduced in length in the presence of disorder, i.e. 
when the system is governed by Hh,e as compared to 
the ordered system described by Hh. However, the x- 
component behaves in a very different manner. De¬ 
pending upon the system parameters, can be both 
higher and lower than its value in the ordered system. In 
Fig.m we exhibit the results in the particular example for 
the system with h/J = 0.3, x = tt/S, which demonstrate 
the random field induced enhancement of mx in the pres¬ 
ence of disorder for two different values of 7 , signaling 
random-field-induced order, also known as “order from 
disorder.” Our numerical observations are further sup¬ 
ported by results obtained analytically via peturvative 
approach at low temperature. The details are discussed 
below. 


A. Perturbative analysis of the low-temperature 
magnetization 


The mean field equations in Eq. ® can be alterna¬ 
tively presented as 



I ar 

Tfi — 

PJ{l + 'y)dmx 

(24) 

and 



(23) 

™ 1 ar 

PJ{1 — 7 ) dmy ’ 

(25) 
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where 

r = Avn [loge Tr exp{-/3Hh)] (26) 

or in the disordered case, 

r = Av^ [logg Tr ex.p{-/3Hh,e)] (27) 

depending on whether the governing Hamiltonian is Hh 
or The symmetry of the Gaussian distribution of 77 

ensures that fh is an even function of e. As a result, 
and vanish at e = 0. Starting with the Eqs. (IMl) 


and (l25]) . straightforward algebra leads to the following 
set of coupled equations: 


drm^ 

de^ 


1 - 


1 


d^r 


Ml+ 7 ) 


d'^edrur 


/3J(l + 7) dml 

d^T d?m, 

H- - 

dmydrrix de^ 


(28) 


and 




1 - 


1 




/3J(1 - 7 ) 


d'^edm,, 


/3J(1 — 7 ) dm'^_ 

j-i 

de^ 


dmydnix 


(29) 


where the total and partial derivatives are evaluated at 

, at e = 0 , 


e = 0. In order to evaluate 77 ^ and 
we need to calculate the partial derivatives at e = 0. For 
example, the expression for ■§^ ^7 is 

(30) 


/3J(1 + 7 ) dm,: 

It follows that 

Av^[(jTOa;(l + 7 ) + ho,){ 


1 


52 ar 


/3J(1 + 7 ) de'^ dmx 
-'iPrf{hy + Jmy{l - 7 ) + eqY 

McMCM) 

rf' tSinh{fik/) 


[cosh (/3fee')] 


/3 


irf‘{hy + Jmy(l — 7 ) + eriY tanh(/3/ce') 


2l3'^rf'{hy + Jmy{l — 7 ) + £ 77 )^ tanh(/3A:£') 
cosh2(/3fc,') 



0.2 0.4 0.6 0.8 I 1.2 
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FIG. 5: (Color online) The span of order-from-disorder in pa¬ 
rameter space. Plots of the functions (a) P{x,j), (b) Q(x,j), 
(c) R{x,j), and (d) S{x,j) with respect to x and j = J/h. 
Note that there are ranges of the (x,j), for which the function 
P(x,j) is positive signaling order from disorder in the system 
described by the Hamiltonian IIh,e- However, this is not true 
for Q(x,j) and R(x,j), which are negative for the entire range 
of X and j. The fact that S{x,j) is also negative in the entire 
range implies that in the presence of disorder, the magneti¬ 
zation vector moves away from the direction of the applied 
random field. All other quantities are dimensionless. Here 
7 = 0.1. 


and 

= (33) 

where the functions P and Q are given by (for j = J/h) 

AE + BC 


and 


Here 


Pix,j) = 


Qix,j) = 


)]• 

(31) 


A{x,j) = acosa:( 


B[x,j) = 6 sina: 


BE - CC 


AC P BD 


DE - CC 


'ib'^ sin^ X 

1 

k'^ 

fc'3 

36^ sin^ X 

3 ' 


7j,/ 5 y^/3 


Here fc/ = + "f)mx) + + ((-1(1 - l)my) + erj ■ 

Next, using the asymptotic expansion of the hyperbolic 
function tanh(,5A:e') Ri 1 — 2 exp(—2/3A:e'), we obtain for 
the partial derivatives at e = 0 

= ( 32 ) 




.. Ja6COScc sinx . 

j) =- ^73 -(1 - T')’ 


Ja& cos X sin x_ , 
C'(x,j) =-p-(I+ 7 ), 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 

























































N >^(1 + 7 )/ oleosa;, , , 

Dix,j) = 1 - ^^(1 - ^ 7 ^), (40) 

and 

£(x,,-)^i-d2^(i-5!|f£), (41) 

with a = 1 + j(l + 7)1 b = 1 + j(l — 7 ), and k' = 
y/a? cos^ X + 6 ^ sin^ x. The positivity of P{x,j) implies 
that order from disorder phenomenon occurs. 

As is clear from Fig. [SKa), there exists a region in the 
parameter space {x,j), for which P{x,j) > 0 which con¬ 
firms that the quenched averaged X-component ma: of 
the magnetization is enhanced by the presence of disor¬ 
der. This does not hold for the quenched averaged Y- 
component, niy, which is reduced in length in the pres¬ 
ence of disorder (see Fig. [SKb)). 

To further investigate the effect of disorder on the 
length m and phase (pi of the magnetization, we expand 
^ as 


TTix \mx ) 


tan<^i = -^1 L.n + (42) 


with 


de^ 


d ruy drrriT. 


l€=0 


le=0 


m% 


le=0 




S{x,j) + 0{e ^), (43) 


where 


S{x,j) = Q{x,j) cosx — P{x,j) sinx. (44) 

S{x,j) is negative for all x and J/h (see Fig. EDd)) im¬ 
plying that the phase always shift towards the X-axis in 
the presence of the random field. 

The square of the length of the magnetization, when 
similarly expanded, is given by 

ml+ml = {ml + m^) | + 2e^ {R{x, j) + 0(e~^)) , 

(45) 

where 


P(x,j) = (Pcosx + Qsmx)j^_^. (46) 

As seen in Fig.EDc), R{x,j) is negative regardless of the 
choice of parameters, i.e., the length of the magnetiza¬ 
tion decreases in the presence of the disorder. Note that 
the analytical results are in agreement with the numer¬ 
ical evidence presented above. It is worth mentioning 
here that the analytical results are valid for small e and 
large pi. The difference between the magnetization in the 
disordered system and the ordered system, as obtained 
analytically, is of the order of e^. Comparison of these 
analytical results with the numerical ones is valid only 
when the same difference, obtained numerically, has pre¬ 
cision of order e^. 


IV. GENERALIZATION TO ARBITRARY 
SPINS AND SCALING OF CRITICAL 
TEMPERATURE 


In this section, our aim is to investigate d-dimensional 
lattices where the occupant of each lattice site is a quan¬ 
tum spin with arbitrary spin angular momentum. Here 
we restrict ourselves to the XX model. For our purposes, 
it is necessary to treat the half-integer and integer spins 
separately. In the following subsections, we derive the 
generalized expressions for the scaling of the magnetiza¬ 
tion and critical temperature for both cases. 


A. Half-integer spins 


The mean-field equations for a general half-integer spin 
(n = 1, 3,...) are 


J-nix Ep=o( 2 p + 1 ) sinh( 2 p -f l)Pk 
^ T,p=o cosh(2p -f l)/3fc 


, (47) 


™ _A \jmy + er]J2p=oi‘^P+^)sinh{2p+l)l3k 

k j:;^,cosH2p+l)pik J’ 

(48) 

where k = yjJ'^ml + {Jmy + eyy)^. 

Finding the magnetization m requires simultaneous so¬ 
lution of the coupled set of Eqs. gZl) and (l4^ . i.e., finding 
the common zeros of the following two functions: 


Fp’^irh) = 

Jmx Ep=o(2p + 1) sinh(2p -f l)j5k 
k Ep=o cosh(2p -f l)j3k 


(49) 


Av^ 


F^y’^im) = 

Jmy + erj Ep=o(2p + 1) sinh(2p -f l)l3k 

k E;=ocosh(2p+l)/3fc J 


(50) 


The Taylor expansion in e, followed by the expansion 
in m, of the functions given in Eqs. (gSl) and (EHl), around 
e = 0 and m = 0, gives 

Fp.’‘‘{fh) = ^[-45 -I- 60J/3s(s + 1) - 

8Js{s + l)(2s^ -I- 2s -I- I)/3^e^]mcos(/>i -I- 

^ gY^/3^ J^[s(s -I- I)(—21(2s^ -I- 2s -I- 1) -I- 

2/3^e^(4s^ -b 2s -I- l)(4s^ -b 6s -b 3 ))]to^ cos(()i -b O(to^), 

(51) 
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FIG. 6: (Color online.) 5p as function of 1/s for the trans¬ 
verse (red circles) and parallel (blue squares) magnetizations 
for e/J = 0.05. The lines serve as guides to the eye. All 
quantities are dimensionless. 


and 

F;’*(m) = 4[-15 -h 20J/3s(s -hi)- 

1D 

8(/s( 5 -|- 1)(25 “h 2s “h j/Tisin 

J^[s(s -h 1)( —21(2s^ -h 2s -h 1) -h 

10/3^e^(4s^ -h 2s -h l)(4s^ -h 6s -h 3 ))]to^ sin^i -h 0(m®), 

(52) 

where s = n-I- 1/2 with n = 0,1,2,.... (j)i has two 
allowed values: 7r/2 (system magnetizes in the direction 
parallel to the disordered field) and 0 (system magnetizes 
in the direction transverse to the disordered field). For 
transverse magnetization, Fy'^(jh) vanishes and two non¬ 
trivial solutions solely come from Eq. as 



[45 - 60J^s(s -h 1) + 8Js(s -h l)(2s2 + 2s + 

J^/3^s{s + l)[-21(2s2 -h 2s + 1) -h 2 / 32 e 2 ( 4 g 2 + 2s -h l)(4s2 -h 6s + 3)]' 


(53) 


The critical point can now be easily obtained by setting 
= 0 in Eq. (l5^ . We get 

8 s(s -h l)(2s2 + 2s + l)jp^e^ - 60s(s + 1) J/3 -h 45 = 0, 

(54) 


which gives 



3 

4Js(s + 1) 


9 (2s2-h2s-hl) 2 
160 J^s^(s -h 1)^ ^ ’ 


(55) 


The critical temperature decreases with the increase of 


number of spin. The shift in critical temperature is of 
the order of for all spins. Note that the generalized 
expressions for the scaling and for the critical tempera¬ 
ture for the pure system with a transverse magnetization 
can be obtained simply by putting e = 0 in Eqs. (I55|) and 
(EH), respectively. 

In order to find the expressions for the parallel magne¬ 
tization, we put (pi = -K12 in the Eqs. (ET|) and (15^ . In 
this case also, the right-hand side of Eq. (EH vanishes to 
leading order, while Eq. (15^ has two nontrivial solutions, 
given by: 


J 


15 - 20Js(s -b l)/3 + 8Js(s -b l)(2s2 -b 2s -b 


8 y J3/33s(s -b l)[-21(2s2 -b 2s -b 1) + 10/32e2(4s2 -b 2s -b l)(4s2 -b 6s -b 3)] ’ 


(56) 


The critical point can be obtained by considering wjj’* = 
0 in Eq. (EHl) and we obtain 


^ 3 27 (2s^ + 2s + 1) 2 

'^’11 4Js(s-b 1) 160 J^s^(s-bl)^ 


(57) 


The generalized expressions of the scaling and the critical 
temperature for the pure system with a parallel magneti¬ 
zation can again be obtained by putting e = 0 in Eqs. (Ei) 
and (E3, respectively. However, the shift in the critical 
temperature due to the random field is in this case big¬ 
ger than in the transverse case and hence the effect of the 


disorder is more prominent in the parallel case similarly 
to what was seen in section (III Bl) without constant field. 


B. Integer spins 


The generalized mean-field equations for the system 
with integer spin n even, are given by 


= Avr, 


Jm, E”=i 2p(e2p/3fc _ e-2p/3fe) 


(58) 
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Jmy + £77 ' 

(59) 

where k = sjJ'^m'^ + {July + erj)^. Now in order to 
find the magnetization m we have to solve the coupled 
set of Eqs. (155)) and OMl). As one can expect from 
the previous discussions, there are two different kinds 
of magnetizations—the transverse magnetization, 
and the parallel magnetization, wjj’*. To derive the crit¬ 
ical scaling for this case, we follow a Taylor expansion 
method, similar to the one used for the half-integer spin 
case. The final expressions for mY‘, the as¬ 

sociated critical temperatures are given by the set of 
Eqs. (1551) . (155]) . (1551) . (1571) with s = n, where = 1, 2, 3,... 

Therefore, we again obtain corrections of order to 
the critical temperature for all the integer spin systems. 
Again, the effect of disorder is more pronounced in the 
parallel magnetization case than in the transverse case. 


rUy = AVy 


C. Critical temperature versus spin quantum 
number 


In order to study the effect of disorder as a function of 
s, we define the dimensionless quantity dp, given by 


< 5/3 



(60) 


dp is shown as a function of 1/s in FiglHlfor e/ J — 0.05. 
We find that the shift in the critical temperature caused 
by the random field decreases with increasing spin quan¬ 
tum number. 


V. CONCLUSIONS 

We considered the quantum spin-1/2 XY model within 
the mean-field approximation and showed that the spon¬ 
taneous magnetization persists in the system with the in¬ 


troduction of a unidirectional quenched disordered field, 
albeit it is smaller than in the pure system. Below a 
certain critical temperature, the magnetization occurs 
in specific directions, either parallel or transverse to 
the disordered field. The critical temperatures and the 
magnitude of the magnetization decrease with increasing 
strength of the disorder. We found perturbative expres¬ 
sions for scaling of the magnetization and the expres¬ 
sions for the scaling of the critical temperatures at which 
the system magnetizes. We also performed numerical 
simulations to obtain the behavior of magnetization for 
various values of the temperature, disorder strength and 
anisotropy parameter, which match with the perturba¬ 
tive calculations for small disorder values. Moreover, we 
extended our analysis to arbitary values of (half-integer 
or integer) spin. We found that the decrease in the length 
of the magnetization due to the random field is of the 
order of the square of the strength of the disorder for 
all values of spin. The magnitude of magnetizations in 
the disordered systems decrease faster with increasing 
spin quatum number and the system requires lower tem¬ 
perature to magnetize when the spin quantum number 
increases, implying that the effect of disorder increases 
with increasing spin quantum number. In addition, we 
studied the random field quantum spin-I/2 XY model 
with an additional constant field, for which we showed a 
random-field-induced ordering in the component of mag¬ 
netization transverse to the disordered field. 
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